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Abstract 

The parentage between Weyl pairs, generalized Pauli group and unitary group is in- 
vestigated in detail. We start from an abstract definition of the Heisenberg-Weyl group on 
the field K and then switch to the discrete Heisenberg-Weyl group or generalized Pauli 
group on a finite ring Z<j. The main characteristics of the latter group, an abstract group 
of order d 3 noted Pd, are given (conjugacy classes and irreducible representation classes 
or equivalently Lie algebra of dimension d? associated with P d ). Leaving the abstract 
sector, a set of Weyl pairs in dimension d is derived from a polar decomposition of SU(2) 
closely connected to angular momentum theory. Then, a realization of the generalized 
Pauli group P d and the construction of generalized Pauli matrices in dimension d are 
revisited in terms of Weyl pairs. Finally, the Lie algebra of the unitary group U(d) is ob- 
tained as a subalgebra of the Lie algebra associated with P^. This leads to a development 
of the Lie algebra of U(d) in a basis consisting of d 2 generalized Pauli matrices. In the 
case where d is a power of a prime integer, the Lie algebra of SU(d) can be decomposed 
into d — 1 Cartan subalgebras. 
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1 Dedicated to the memory of my teacher and friend Moshe Flato on the occasion of the tenth anniversary 
of his death. 



1 Introduction 

The present paper is devoted to three major ingredients of quantum mechanics, namely, 
the Heisenberg-Weyl group connected with Heisenberg commutation relations [1], the 
Pauli spin matrices [2] used in generalized angular momentum theory and theory of uni- 
tary groups, and the pairs of Weyl [3] of relevance in finite quantum mechanics. 

The Heisenberg-Weyl (or Weyl-Heisenberg or Heisenberg) group HWQ3L), also called 
the Weyl group 01, is of central importance for the quantization process and its Lie al- 
gebra turns out to be a basic building unit for quantum mechanics [5|. Note that the 
Lie algebra of HW(M) should not be confused with the Weyl-Heisenberg algebra (or 
oscillator algebra spanned by the creation, annihilation and number operators) and its 
supersymmetric extensions Wk [6]. 

A discrete restriction HW(Zd) of HW(R), corresponding to the replacement of the 
infinite field R by a finite ring Zj = Z/dZ, yields a group of order d 3 (d arbitrary in 
N \ {0, 1}). This group was introduced by Stovicek and Tolar [7] in connection with 
quantum mechanics in a discrete space-time, by Balian and Itzykson in connection with 
finite quantum mechanics [8], and by Patera and Zassenhaus [9] in connection with grad- 
ings of simple Lie algebras of type A n _ 1 . The case where the ring Z d is replaced by a 
finite (Galois) field ¥ g gave rise to several mathematical studies ifTOl ITTTl . The discrete 
Heisenberg-Weyl group, also known as the generalized Pauli group, plays a central role 
in quantum information, cf. the interest of Galois fields in finite quantum mechanics lfT2l 
and, consequently, in quantum information and quantum computation. In this connection, 
a finite Heisenberg-Weyl group was used for a description of phase oscillations of EPR 
states IPT3H . 

What is the relationship between the Heisenberg-Weyl group and Weyl pairs? First 
of all, a definition of a Weyl pair is in order. A Weyl pair (X, Z) in d dimensions is a 
pair of rf-dimensional unitary matrices X and Z that satisfy the g-commutation relation 
XZ — qZX = and the cyclic relations X d = Z d = I (I standing here for the unitary 
matrix), where q is a primitive root of unity with q d = 1. The concept of a pair of Weyl, 
initially introduced for dealing with quantum dynamical systems in finite dimension ||3), 
was used for the construction of unitary bases in finite-dimensional Hilbert spaces [fT4l 
and (independently) for the factorization of the secular equation corresponding to finite- 
dimensional eigenvalue problems [fT5l . In the last 20 years, the notion of Weyl pairs was 
used for the construction of generalized Pauli matrices in domains as different as graded 
Lie algebras and quantum information. 

The usual Pauli matrices a x , cr y and a z are useful for the representation theory of the 
Lie group SU(2). Therefore, a natural extension of the Pauli matrices resulted in the 



sixties from the interest of the group 577(3) for the classification of elementary particles 
IfToll . This gave gise to the Gell-Mann matrices and the Okubo matrices. Further exten- 
sions of the Pauli matrices came out of the introduction of the group S77(4) for charmed 
particle IfTTl and of the group S77(5) for a grand unified theory of quarks and leptons 
|[T8l . The Gell-Mann lambda matrices for S77(3) and their extension to Cartan bases for 
SU(d) undoubtedly constitute a systematic extension of the ordinary Pauli matrices. This 
statment is particularly justified as far as the tensor structure (involving symmetric and 
antisymmetric tensors) of their algebra is concerned [19]. We shall deal in this paper with 
another extension of the Pauli matrices in d dimensions which turns out to be of special 
interest in the case where d is a power of a prime integer. Indeed, generalized Pauli matri- 
ces can be constructed in a systematic way by making use of Weyl pairs. In this direction 
let us mention the pioneer work of Patera and Zassenhaus [9[. In the last two decades, 
the construction of generalized Pauli spin matrices has been extensively used in the the- 
ory of semi- simple Lie algebras, in quantum mechanics (complete state determination, 
reconstruction of a density matrix, discrete Wigner functions), in quantum information 
and quantum computation (mutually unbiased bases, unitary error bases, quantum error 
correction, random unitary channels, mean king's problem, positive operator valued mea- 
sures, and quantum entanglement), and in the study of modified Bessel functions (see for 
instance [HI! and [|20ll-|]32l). 

From a group-theoretical point of view, the rf-dimensional generalized Pauli matrices 
may serve to construct a generalized Pauli group in d dimensions, a group generalizing 
the ordinary Pauli group spanned by the ordinary Pauli matrices (see lT71l- lfT2"l . ETTl . [|24ll 
and [1331 - [1421 ). In fact, this group is nothing but the discrete Heisenberg-Weyl group 
HW(Zd). This generalized Pauli group has been recently the object of numerous studies 
partly in connection with the Clifford or Jacobi group 11341 [361 [371 [381 as well as graph- 
theoretical and finite-geometrical analyses of the generalized Pauli operators 14111421 . 

The object of this work is to further study the link between the Heisenberg-Weyl 
group, the Weyl pairs, the generalized Pauli matrices and the generalized Pauli group and 
to revisit their interest for unitary groups. We shall start with an abstract definition of 
the Heisenberg-Weyl group, pass to an abstract version of HW(Zd) an d briefly study it. 
Then, we shall deal with the introduction of Weyl pairs from a polar decomposition of the 
Lie algebra su{2) and we shall use them for finding a realization of HW(Z d ) isomorphic 
to the generalized Pauli group in d dimensions. Finally, some of the generators of the Pauli 
group in d dimensions shall be used for constructing the Lie algebra su(d) of SU(d) in a 
basis that is especially adapted, when d is a power of a prime integer, to a decomposition 
of su(d) into a direct sum of d + 1 Cartan subalgebras. 



2 The Heisenberg-Weyl group 
2.1 The Lie group HW (I 



We start with an abstract definition of the Heisenberg-Weyl group HW(M). Let us con- 
sider the set of triplets 

S := {(x,y, z) : x,y, z E R}. (1) 

The set S can be equipped with the internal composition law S x S — ► S defined trough 

(x,y,z)(x',y',z') := (x + x' - zy',y + y',z + z'). (2) 

It is clear that the set S is a group with respect to the law ©. We denote HW(M) this 

group and call it the Heisenberg-Weyl group (for evident reasons to be given below) on 

the infinite field E. More precisely, we have the following result. 

Proposition 1. The group HW(R) is a noncompact Lie group of order 3. This nonabelian 

group is nilpotent (hence solvable) with a nilpotency class equal to 2. 

Proof. The proof is trivial. Let us simply mention that the nilpotency of HW(K) follows 

by repeated use of the commutator 

[<y, y\ *')» (x, y, z)\ = (zy' - yz', 0, 0) (3) 

of the elements (x f , y', z') and (x, y, z) of the group HW(R). Equation © shows that 
(x,y,z) and (x', y',z') commute if and only if zy' — yz 1 = 0. □ 

In the terminology of Wigner 11431 , the group HW(R) is not ambivalent (ambivalent 
means that each conjugacy class contains its inverse elements). Indeed, since 

(x,y, z)' 1 = (-x-yz, -y, -z) (4) 

and 

(x',y',z')(x,y,z)(x',y',z')~ 1 = (x + zy' - yz',y,z) (5) 

it is evident that only the class C( ,o,o) = {(0, 0, 0)} of the identity element (0, 0, 0) is 
ambivalent. 

2.2 The Lie algebra of HW(R) 

We may ask why to call HW(M) the Heisenberg-Weyl group? The following result clar- 
ifies this point. 



Proposition 2. A set of infinitesimal generators of HW(M) is 

19 1 d _ 1 f d d\ 

i dx i dy i \dz dx J 

This set of generators satisfies the formal commutation relations 

[Q,P]_=iH [P,H]- = [H,Q]. = Q (?) 

with H = H, Q = Q and P = V. The Lie algebra hw(R) of HW(R), with the Lie 
brackets ((7]), is a three-dimensional nilpotent (hence solvable) Lie algebra with nilpotency 
class 2. 

Proof. The proof easily follows by working in a neighbourhood of the identity (0, 0, 0) of 
HW(M) and by considering the series w 1 = hw(M), w 2 = [w 1 , u* 1 ]., w 3 = [w 1 , w 2 ]_, ■ • • 
where [A, B]_ refers here to the set of commutators [«,/?]_ with a e A and /3 e B. □ 

The connection with the Heisenberg commutation relations is clearly emphasized 
by ©. This constitutes a partial justification for calling HWQSL) the Heisenberg-Weyl 
group on IR. The Lie algebra hw(M) was derived from a matrix group [4] and stud- 
ied at length from the point of view of quantum mechanics Q. This algebra admits 
infinite-dimensional representations by Hermitean matrices. In particular, we have the 
infinite-dimensional harmonic oscillator representation which is associated with the op- 
erator realization H = Hho '■= hi, Q = Qho '■= x and P = Vho '■= f ■§-, where h is the 
rationalized Planck constant. On the other side, we may expect to have finite-dimensional 
representations of hw(R) at the price to abandon the Hermitean character of the represen- 
tation matrices. 

As an example, we have the three-dimensional representation of hw(R) defined by 
H = H 3 ,Q = Q 3 and P = P 3 with 
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H 3 := Q 3 := \i P 3 := . (8) 

\i 0/ \0 0/ \0 -i 0/ 

We can look for the matrix Lie group which corresponds to the Lie algebra spanned by 
the set {H 3 , Q 3 , P 3 }. This yields Proposition 3. 
Proposition 3. The exponentiation 

M(x, y, z) := exp[i{xH 3 + yQ 3 + zP 3 )} (9) 

leads to 



M(x,y,z)= -y 10. (10) 




The matrices M(x, y, z) satisfy the composition law 

M(x, y, z)M(x', y', z) = M(x + x + -zy' - -yz\ y + y',z + z) (11) 

so that the set 5" := {M(x, y, z) : x,y, z e IR} endowed with the law (fTTj) is a group 
isomorphic to HW(R). 

Proof. A simple expansion of © where H 3 , Q 3 and P 3 are given by © yields (flOl) . The 
isomorphism follows from the bijection S — > S" : (x,y, z) i— > M(—x — \yz, —y, —z). 
Note that the matrix form (flOl) corresponds to two other sets {7^±, Q±, P±} of infinitesi- 
mal generators of HW(M), namely, 

dx ~"~ \dy 2 dx ) ' \dz 2 dxj 

which satisfies © with H = H±, Q = Q± and P = V± (cf. 0H3). □ 



3 The Pauli group 

3.1 The abstract Pauli group 
3.1.1 The group P d 

We shall be concerned in this section with a discretization of the Heisenberg-Weyl group 
HW(J3L). A trivial discretization of HW(J3L) can be obtained by replacing the field IR 
by the infinite ring Z. This leads to an infinite-dimensional discrete group HW(Z). A 
further possibility is to replace IR by the finite ring Z^ = Z/dZ where d is arbitrary in 
N \ {0, 1}. (In the case where d is a prime p or a power of a prime p e with e G N \ {0, 1}, 
the finite ring Z/dZ can be replaced by the Galois field F p or ¥ p e .) This yields a finite 
group HW(Z d ) which can be described by the following result. 
Proposition 4. The set 

S d := {(a,b, c) : a, b, c e Z d } (13) 

with the internal composition law S d x S d —> S d defined trough 

(a, b, c){a', b', c') := (a + a! - cb\ b + b',c + d) (14) 

(where from now on the addition is understood modulo d) is a finite group of order d 3 . 
This nonabelian group HW(Z d ), noted P d for short, is nilpotent (hence solvable) with a 
nilpotency class equal to 2. 



Proof. The proof of Proposition 4 is elementary. Note simply that we have the canonical 
decomposition 

(a, b, c) = (a, 0, 0)(0, 6, 0)(0, 0, c) (15) 

for any element (a, 6, c) of P d and that two elements (a, b, c) and (a', b', d) of P,i commute 
if and only if cb' — be' = (mod d). □ 

We call the abstract group P d the (generalized) Pauli group in d dimensions. At this 
stage, we can give the main reason for associating Heisenberg, Pauli and Weyl in the title 
of the present paper. As a point of fact, the discretization of the group HW(M), a group 
associated with the Heisenberg commutation relations, via the replacement R — > Z/dZ 
gives rise to the group P d , a group which can be realized in terms of generalized Pauli 
matrices, which in turn can be constructed in terms of Weyl pairs (see below). 

3.1.2 Some subgroups of P d 

Among the subgroups of Pd, we can mention proper subgroups of order d and d 2 (there 
are no other proper subgroups if d is a prime integer). We simply list below the subsets of 
Sd, which together with the law (fT4l) . provide us with some important subgroups of Pd. 

- The set {(a, 0, 0) : a e Z d } gives an invariant abelian subgroup of Pd of order d 
isomorphic to the cyclic group 7h d . In fact, this subgroup is the centrum Z(P d ) of P d and 
P d /Z(P d ) is isomorphic to Z rf ® Z rf . 

- The set {(0, b, 0) : b e Z d } gives an abelian subgroup of Pd of order d isomorphic to 
Z d . 

- Similarly, the set {(0, 0, c) : c G Z d } gives also an abelian subgroup of P d of order d 
isomorphic to Z^. 

- The sets {(a, b, 0) : a,b E Z d } and {(a, 0, c) : a,c <E Z, d } give two invariant abelian 
subgroups of P d of order c? 2 isomorphic to Z d ® Z d . 

- Finally, the set {(a, 6, b) : a, 6 e Z d } give an invariant abelian subgroup of P d of 
order d 2 . 

3.1.3 Conjugacy classes of P d 

The conjugacy classes of P d readily follow from 

(a , 6', c')(a, 6, c)(a', b\ c') _1 = (a + c6' — be , b, c) (16) 

with addition modulo d. This can be precised by the following result. 

Proposition 5. The group P d has d(d + 1) — 1 conjugacy classes: d classes containing 

each 1 element and d 2 — 1 classes containing each d elements. 
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Proof. It can be checked that the class C( a ,o,o) °f ( a , Oj 0) is £(a,o,o) = {(a, 0, 0)}; therefore, 
there are d classes with 1 element. Furthermore, the class C( a ,&, c ) °f ( a ; &j c )> w i m m e case 
b = c = excluded, is C( ai &, c ) = {(a', 6, c) : a' G Z rf }; this yields d 2 — 1 classes with d 
elements. We note that the group Pd is not ambivalent in general. □ 

The case d = 2 is very special since the group P 2 of order 8 is ambivalent like the 
group Q of ordinary quaternions, another group of order 8. Not all the subgroups of P 2 are 
invariant. Therefore, the group P 2 is not isomorphic to Q (for which all subgroups are in- 
variant). Indeed, it can be proved that P 2 is isomorphic to the group of hyperbolic quater- 
nions associated with the Cayley-Dickson algebra A(ci,c 2 ) with (ci,c 2 ) ^ (—1,-1) 
defined in Il44ll . In this respect, the Pauli group P 2 defined in this work differs from the 
Pauli group in d — 2 dimensions considered by some authors, a group isomorphic to the 
group Q of ordinary quaternions. Let P 2 be this latter Pauli group. It consists of the 
elements a : = ±cr , ±io x , ±ia y , ±ia z (where a is the 2x2 unit matrix). Let us also 
mention that an extension of the group P 2 is used in quantum computation P31 (see also 
If33ll42l ). This extension, say P 2 , is obtained from a doubling process: The group P 2 
consists of the elements of the set {a, ia : a G P 2 }. Thus, the conjugation classes and the 
irreducible representation classes of P 2 trivially follow from those of P' 2 . 

3.1.4 Irreducible representations of Pd 

The duality between conjugacy classes and classses of irreducible representations leads 

to the following result. 

Proposition 6. The group Pd has d(d + 1) — 1 classes of irreducible representations: d 2 

classes of dimension 1 and d— 1 classes of dimension d. 

Proof. It is sufficient to apply the Burnside-Wedderburn theorem. □ 

As a corollary of Proposition 5 and Proposition 6, the difference between the order 
of Pd and its number of classes (conjugacy classes or irreducible representation classes) 
is odd if d = 2k (k G W), or a multiple of 16 if d = Ak + 3 (k G N) or a multiple of 
32 if d — Ak + 1 (k G N*). (For an arbitrary finite group of odd order, the difference 
is a multiple of 16.) Furthermore, the number of elements of Pd which commute with a 
given element (a, b, c) of Pd is d 3 or a multiple of d 2 according to whether the order of 
the conjugation class containing (a, b, c) is 1 or d; see [HTI for a more elaborated result, in 
the form of a universal formula, and its interpretation in terms of the fine structure of the 
projective line defined over the modular ring Z d . Note that Proposition 5 and Proposition 
6 are in agreement with the results obtained [10] in the case where d is a power of a prime 
integer corresponding to the replacement of the ring Z d by the Galois field ¥ d . 



3.1.5 A Lie algebra associated with P d 

We close the study of the abstract group Pd with a result devoted to the association of 
Pd with a Lie algebra of dimension d 3 . Let us consider the group algebra (or Frobenius 
algebra) F(Pd) of the generalized Pauli group Pd. Such an algebra is an associative 
algebra over the field C. By applying the process developed by Gamba [|46l . we can 
construct from F{Pd) a Lie algebra, which we shall denote as pd, by taking 

((a, b, c), (a', b' } c')) := (a + a' - cb', b + b',c + c') - (a + a' - be', b + b',c + c') (17) 

for the Lie bracket of (a, b, c) and (a', b', c'). (The right-hand side of (flTT) is defined in 

F(P d ).) The set Sd constitutes a basis both for the Frobenius algebra F(P d ) and the Lie 

algebra p d (S d is a Chevalley basis for p d ). As a further result, we have the following 

proposition. 

Proposition 7. The Lie algebra pd of dimension d 3 , associated with the finite group P d of 

order d 3 , is not-semi-simple. It can be decomposed as the direct sum 

d? d-l 

Pd = M (l)0 M (d) (18) 

i i 

which contains d 2 Lie algebras isomorphic to u(l) and d — 1 Lie algebras isomorphic to 
u(d). 

Proof. The proof can be achieved by passing from the Chevalley basis of pd, inherent to 
(fT7l) . to the basis generated by the idempotent (or projection) operators and nilpotent (or 
ladder) operators, defined in F(Pd), associated with the classes of irreducible representa- 
tions of Pd- Equation ([TBI is reminiscent of the fact that Pd has d 2 irreducible representa- 
tion classes of dimension 1 and d — 1 irreducible representation classes of dimension d. 
D 

3.2 A realization of the Pauli group 
3.2.1 Polar decomposition of SU(2) 

Let £(2j + 1), with 2j E N, be a (2j + 1) -dimensional Hilbert space of constant angular 
momentum j. Such a space is spanned by the set { \j, m) : m = j, j — 1, • • • , — j}, where 
\j, m) is an eigenstate of the square j 2 and the z-component j z of a generalized angular 
momentum iPJTl . The state vectors \j, m) are taken in an orfhonormalized form, i.e., the 
inner product (j, m\j', m') is equal to 5 m , m '. 

Following the approach of 11481 . we define the linear operator v ra via 

v ra \j, m) = (1 - S md ) g^- m)a \j, m + 1) + 6 mj e i2 ^ r \j, -j) (19) 



where 



r G 



a e 



q = exp 



2vri 



& + 1 
The matrix V ra of the operator v ra in the spherical basis 

reads 



(20) 



(21) 



I',, 



/ <f 

q 2a 



V e i2wjr 



\ 



J2ja 

1 

o I 



(22) 



The matrix V ra constitutes a generalization of the matrix V a introduced in [|49l (see also 
EQl). 

The shift operator v ra takes its origin in the study of the Lie algebra of 577(2) in a 
nonstandard basis with the help of two quon algebras describing g-deformed oscillators 
[50J. The operator v ra is unitary. Furthermore, it is cyclic in the sense that 



v2j+l _ e i2TTJ(a+r)j 



(23) 



where I is the identity. The eigenvalues and eigenvectors of v ra are given by the following 

result. 

Proposition 8. The spectrum of the operator v ra is nondegenerate. For fixed j, r and a, it 

follows from 



v ra \ja;ra) = q^ a+ ^ a \ja;ra) 



(24) 



where 



\ja; ra) 



VW+ 



1 j 

— y 



(j+m)(j—m+l)a/2—jmr-\-(j+m)a 



\j, m) 



(25) 



m=-j 



for a = 0, 1, • • • ,2j. 

A second linear operator is necessary to define a polar decomposition of S77(2). Let 
us introduce the Hermitean operator h through 



h\j,m) = y/(j + m)(j -m + l)\j,m). 



(26) 
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Then, it is a simple matter of calculation to show that the three operators 

j+ = hv ra 3- = vl a h j z = -(h 2 - vl a h 2 v ra ) (27) 

satisfy the ladder equations 

j+\j,m) = g +(j - m+s - 1/2) V(j -m)(j+m + l)|j, m + 1) (28) 

j-\j,m) = g-°'- m+s+1/2) V(i + m)(j-m + 1) \j, m - 1) (29) 

and the eigenvalue equation 

j z \j,m) =m\j,m) (30) 

where s = 1/2. (Note that there is one misprint in the corresponding relations of [|40l .) 

Therefore, we have the following result. 

Proposition 9. The operators j + , j_ and j z satisfy the commutation relations 

[jz ,j+] = +J+ [Jz,J-] = -J- \j+ , 3- } = 2 Jz (3 1 ) 

and thus span the Lie algebra of SU(2) over the complex field. 

The latter result does not depend on the parameters r and a. However, the action of 
j + and j_ on \j,m) on the space £(2j + 1) depends on a (an a priori real parameter 
to be restricted to integer values in what follows); the usual Condon and Shortley phase 
convention used in spectroscopy corresponds to a = 0. The writing of the ladder operators 
j + and j_ in terms of h and v ra constitutes a two-parameter polar decomposition of the 
Lie algebra of SU(1, 1) [or SU (2)over the complex field]. This decomposition is an 
alternative to the polar decompositions obtained independently in ll5Tll52l . 

3.2.2 Weyl pairs 

The linear operator x := v 00 such that (cf. (fT9l) ) 

x\j, m) = (1 - S m ,j) \j, m + 1) + 6 m j\j, -j) (32) 

has the spectrum (1, q, ■ • • , g 2jf ) on £(2j + 1). Therefore, the matrix X := Voo of x on the 
basis b s is unitarily equivalent to 

Z:=diag(l,g,---,g 2j ). (33) 

The linear operator z corresponding to the matrix Z can be defined by 

z\j,m)=qi- m \j,m). (34) 

11 



The two isospectral operators x (a cyclic shift operator) and z (a cyclic phase operator) 
are unitary and constitute a pair of Weyl (x, z) since they obey the g-commutation relation 

xz — qzx = (35) 

(or XZ — qZX = in matrix form). These two operators are connected via 

x = pzf & z = fxf (36) 

where / is the Fourier operator such that 

f\j,m) = —±— Y, q- {j - m){j - m,) \j,m'). (37) 

m'=—j 

The operator / is unitary and satisfies 

/ 4 = 1 (38) 

(see d52l for a general treatment of Fourier operators on finite-dimensional Hilbert spaces). 
Let F be the matrix of the linear operator / in the basis b s . Indeed, F is a circulant ma- 
trix. Note that the reduction of the endomorphism associated with the matrix X yields the 
matrix Z. In other words, the diagonalization of X can be achieved with the help of the 
matrix F via Z = FXFl 

We conclude that the polar decomposition of SU(2) described in Section 3.2.1 pro- 
vides us with an alternative derivation of the Weyl pair (X, Z) . Of course, other pairs 
of Weyl (V ra , Z), corresponding to (v ra , z) with the property v ra z — qzv ra = 0, can be 
derived by replacing v 00 by v ra . Note that v ra = v r0 z a . 

3.2.3 Weyl pairs and Pauli group 

Let us define the d 3 operators 

w abc := q a x b z c a,b } ce 1 d . (39) 

The action of w abc on the Hilbert space £ (2j + 1) is described by 

w abc \j,m)=q a+ ^ c \j,m + b) (40) 

where m + b is understood modulo 2j + 1. The operators w a bc are unitary and satisfy 

Tr£-( 2 j+i) (wl bc w a , b > c ,) = q a '~ a d 5 bjbf 5 C<C > (41) 
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with d := 2j + 1. In addition, we have the following central result. 

Proposition 10. The set Wd := {ww : a,b,c E Z d } endowed with the multiplication of 

operators is a group isomorphic to the Pauli group P,i. Thus, the group Pd is isomorphic 

to a subgroup of U(d) for d even or SU(d) for d odd. 

Proof. The proof is immediate: It is sufficient to consider the bijection Wd — ► Sd '■ 

Wa.bc l— > (°) &> c )> to use repeatedly (|35l) or (|40l) . and to note that the matrix of w a fe c in the 

basis b s belongs to U(d) for d even and to SU(d) for c/ odd. As a consequence, the Lie 

bracket ((a,b,c), (a',b',c')), see (fTTT) . corresponds to the commutator [iu a & c , uvftv]- so 

that the Lie algebra pd associated with the finite group Pd corresponds to the commutation 

relations 

[Wabc, «V&/ c /]_ = W a(3l - W a >p>y (42) 

with a = a + a' — cb', (3 = b + 6', 7 = c + d , a' = a + cb' — be', j3' = (3 and 7' = 7. □ 

3.2.4 Weyl pairs and infinite-dimensional Lie algebra 

We close this section by mentioning another interest of Weyl pairs (v ra , z). By defining 
the operators 

t m = q \^rn 2v rrn z m 2 m = (jn^ m 2 ) G N* 2 (43) 

we easily obtain the following result. 

Proposition 11. The commutator of the operators t m and t n reads 



[t m ,t n }_ = 2isin I ', 



-m A n t m+n (44) 

VA7 t- J 

where 

m An = vf\\U2 — 77^2^1 m + n = [m\ + m, m2 + n%). (45) 

Therefore, the linear operators t m span an infinite-dimensional Lie algebra. 

The so-obtained Lie algebra is isomorphic to the algebra introduced in ll53~l . The latter 
result parallels the ones derived, on the one hand, from a study of A;-fermions and of the 
Dirac quantum phase operator through a g-deformation of the harmonic oscillator ll54l 
and, on the other hand, from an investigation of correlation measure for finite quantum 
systems [|55l . 
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3.3 Mutually unbiased bases 

We now briefly establish contact with quantum information. For this purpose, let us 
introduce the notation 

k-=j-m \k):=\j,m) d:=2j + l. (46) 

Thus, the angular momentum basis {\j, j), \j, j — 1), • • • , \j, —j)} of the finite-dimensional 
Hilbert space £(2j + 1) reads {|0), |1), • • •, \d— 1)}. Let us note 

B d :={|A;>:A;=0,l,"-,d-l} (47) 

the latter orthonormal basis, known as the computational basis in quantum information 
and quantum computation. From now on, the real number a occurring in (|25b shall be 
restricted to take the values a = 0,1, ■ ■ ■ ,d— 1. 

From equation (|25l ), we can write the eigenvectors \aa) := \ja; Qa) of the operator 

w 0a as 

d-i 

Vd 



\aa) = * ^ g (d-*-i)(*+i)«/2-(*+i)«|fc) (48) 



fe=0 

where, for fixed a (a — 0, 1, • • • , d — 1), the index a takes the values 0, 1, • • • , d — 1. Note 
that 

So«:={|oa):a = 0,l,---,d-l} (49) 

is another orthonormal basis of £ (d). 

Proposition 8 can be transcribed in matrix form by using the generators E XjV of GL(d, C) 
(see also [40] where a different normalization is used). The d x d matrix E xy (with 
x, y G Zrf) is defined by its matrix elements 

(E X} y) M = S k , x Si, y k, I G Z d . (50) 

Therefore, the matrix V 0a of the operator v 0a in the computational basis B d is 

d-2 



V 0a = E d . lfi + J2 q {k+1)a E k , k+ x. (51) 

fc=0 

The eigenvectors <p{aa) of the matrix V 0a are expressible in terms of the d x 1 column 
vectors e x (with x G Z d ) defined via 

( e *)fco = ^,^' k e Z d- (52) 
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In fact, we can check that 

^aa) = ^~y g (*-*-i)(*+i)«/2-(fc+i)« efc (53) 

^ d t^o 

satisfies the eigenvalue equation 

V 0a <p(aa) = q {d - l)a/2 - a ^{aa). (54) 

Furthermore, the d x d matrix 

H " '■= E E g( d - fe - 1 )( fc + 1 ) a / 2 -( fe + 1 ) a J B fejQ (55) 

«=0 fc=0 

reduces the endomorphism associated with V 0a . In other words, we have 

d-l 

H% a H a = q^ 1 ^ 2 d Y, q' a E a , a . (56) 

Note that H a is a generalized Hadamard matrix in the sense that 

H\H a = dl (57) 

and the modulus of any element of H a is unity. Observe that the Fourier matrix F can be 
written as 

1 d-X 

F = (H Sy S:=^=J2E Ptd _ (58) 

Vd^ Q 

where S acts as a permutation matrix normalized by 4=. 

As an application of (l48l or (l53l) to mutually unbiased bases, we have the following 
result (see also ||40l|49l). 

Proposition 12. In the case where d = p is a prime integer, the bases B 0a for a = 
0, 1, ■ ■ • , p — 1 together with the computational basis B d constitute a complete set of p + 1 
mutually unbiased bases. 
Proof. According to the definition of mutually unbiased bases 11561 . we need to prove that 

\{k\aa)\ = ^- (59) 



and 



\{aa\b(3)\ = 5 a> p5 a , b + — (1 - 5 , 6 ) (60) 
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for any value of a, b, a, (3 and k in Z^. Equation (1591 ) simply follows from (1481) and 
equation (|60l ) was proved in [|40ll by making use of generalized quadratic Gauss sums. □ 

The interest of (|48l) or (|53l) with d = p,p prime (including the case p = 2), is that the 
p 2 vectors corresponding to the p mutually unbiased bases besides the computational basis 
are obtainable from one single formula that is easily codable on a computer (the single 
formula corresponds to the diagonalization of only one matrix, namely, the matrix Vo a 
where a can take the values a = 0, 1, • • • , p — 1). In matrix form, the p mutually unbiased 
bases besides the computational basis are given by the columns of the Hadamard matrices 
matrices H a (a = 0, 1, • • • ,p — 1). 

Going back to d arbitrary, we can check that the bases B 00 , B m and B d constitute a 
set of 3 mutually unbiased bases. Therefore, we recover a well-known result according to 
which there exists a minimum of 3 mutually unbiased bases when d is not a prime power. 

4 Weyl pairs and unitary group 

In this section, we shall focus our attention on one of the u(d) subalgebras of p d . Such 
a subalgebra can be constructed from a remarkable subset of {ww : a,b,c £ Z d }. This 
subset is made of generalized Pauli operators. It is generated by the Weyl pair (x, z) or 
(X, Z) in matrix form. 

4.1 Generalized Pauli operators 

Following the work by Patera and Zassenhaus 0, let us define the operators 

Uab ■ = w ab = x a z b a, b e Zd- (61) 

The operators u a b are unitary. Note that the matrices X a Z h of the operators u a t, in the basis 
b s belong to the unitary group U(d) for d even or to the special unitary group SU(d) for d 
odd. The d 2 operators u a b shall be refered to as generalized Pauli operators in dimension 
d. It should be mentioned that matrices corresponding to the operators of type (|6TT) were 
first introduced long time ago by Sylvester l(57l in order to solve the matrix equation 
PX = XQ; in addition, such matrices were used by Morris ll58l to define generalized 
Clifford algebras in connection with quaternion algebras and division rings. The operators 
u a b satisfy the two following properties which are direct consequences of (|4TT) and (1421) . 
Proposition 13. The set {u a b : a,b £ Z d } is an orthogonal set with respect to the Hilbert- 
Schmidt inner product. More precisely 

Tr£(2i+1) (UabUa'b'J = d 5 a ,a> <W (62) 
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where the trace has to be taken on the d-dimensional space £(2j + 1) with d :— 2j + 1. 
Proposition 14. The commutator [u ab , u a ty\- and the anti-commutator [u ab , u a iy] + ofu ab 
and u a iy are given by 

[u ab , u a , b ,} T = (q~ ba ' T q- aV ) u a , lb „ a" := a + a' b" := b + b'. (63) 

Consequently, [u ab , Ua/v]- = if and only if ab' — ba! = (mod d) and [u abl u a 'v]+ — 
if and only if ab' — ba' = (l/2)d (mod d). Therefore, all anti-commutators [u a b, u a 'b']+ 
are different from if d is an odd integer. 

The d 2 pairwise orthogonal operators u ab can be used as a basis of the Hilbert space 
C d (with the Hilbert-Schmidt scalar product) of the operators acting on the Hilbert space 
C d (with the usual scalar product). In matrix form, they give generalized Pauli matrices 
in (2j + 1) x (2j + 1) dimensions, the spin angular momentum j = 1/2 corresponding 
to the ordinary Pauli matrices. 

Example 1: j = 1/2 =>■ q = — 1 and d — 2. The matrices of the four operators u ab with 
a, b = 0, 1 are 

'0 <70 _ / 1 u \ _ 1 



X^ u =( J) X = X^=K ;) (64) 



-0 7 l /I \ v . v-1 71 






n 


1 


oj 





-1 


1 






^ = •^=^0 -ij 1':=.^'=^ j. (65) 

In terms of the usual (Hermitean and unitary) Pauli matrices a x , a y and a z , we have X = 
a x ,Y = —ia y , Z = a z . Note that a normalization for the Pauli matrices different from the 
conventional one is also used in [0. The group-theoretical approaches developed in [|9l 
and in the present paper lead to Pauli matrices in dimension 2x2 that differ from the usual 
Pauli matrices. This is the price one has to pay in order to get a systematic generalization 
of Pauli matrices in arbitrary dimension (see also [|9l |23l). It should be observed that 
the commutation and anti-commutation relations given by (|63l with d — 2 correspond to 
the well-known commutation and anti-commutation relations for the usual Pauli matrices 
(transcribed in the normalization X 1 Z° = a x , X 1 Z 1 = —ia y , X°Z 1 = a z ). 
Example 2: j = 1 =^> q = exp(27ri/3) and d = 3. The matrices of the nine operators u a b 
with a, b — 0, 1,2, viz., 

X°Z° = I X Y Z° = X X 2 Z° = X 2 X°Z 1 = Z X°Z 2 = Z 2 (66) 
X X Z X = XZ X 2 Z 2 X 2 Z X = X 2 Z X X Z 2 = XZ 2 (67) 
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are 





X = 1 X 2 = 1 (68) 






Z 2 = I q 2 I XZ =- I g 2 I (69) 






X 2 Z 2 =1 X 2 Z =1 XZ 2 =0 g . (70) 




These matrices differ from the Gell-Mann matrices lfT6ll used in elementary particle physics. 
They constitute a natural extension of the Pauli matrices in dimension 3x3 (see also 

4.2 The unitary group in the generalized Pauli basis 

From Proposition 14, it is clear that the set {u ab : a,b = 0, 1, • • • , d — 1} can be used 
as a set of generators of the Lie group U(d). Thus the generalized Pauli matrices X and 
Z form an integrity basis for the Lie algebra of U(d). This can be precised by the two 
propositions below. 

Proposition 15. The set {X a Z b : a, b — 0, 1, • • • , d — 1} form a basis for the Lie algebra 
u(d) of the unitary group U(d) for d arbitrary. The Lie brackets of u(d) in such a basis 
(that we denote as the Pauli basis) are given by 

[X a Z\ X a 'Z b '\_ = J2( ab > a ' h 'i a"b")X a "Z b " (71) 

a" b" 

where the structure constants (ab, a'b'; a"b") read 

(ab, a'b'; a"b") = 6{a", a + a')5(b", b + V) (q- ba ' - q- ab '^ (72) 

with a,b,a',b' = 0, 1, • • • , d — 1 (mod d). The structure constants (ab, a'b'; a"b") with 
a" = a + a' and b" = b + V are cyclotomic polynomials associated with d. They vanish 
for ab' — ba! = (mod d). 

Proposition 16. In the case where d = p is a prime integer, the Lie algebra su(p) of 
the special unitary group SU(p) can be decomposed into a direct sum of p + 1 abelian 
subalgebras of dimension p — 1. More precisely 

su{p) =v Q ®v 1 ®---®v p (73) 
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where each of the p + 1 subalgebras v , t>i, • ■ • , v p is a Cartan subalgebra generated by a 
set of p — 1 commuting matrices. The various sets are 



V, 



Vi 


= {x 1 ^ ,^ 2 ^ ,---,^- 1 ^ } 


(74) 


v 2 


= {x l z\x 2 z 2 ,---,x p - 1 z p - 1 } 


(75) 


v 3 


= {x 1 z 2 ,x 2 z\---,x p - 1 z p - 2 } 


(76) 
(77) 


>-l 


= {x l z p - 2 ,x 2 z p -\---,x p - 1 z 2 } 


(78) 


v P 


= {x l z p -\x 2 z p - 2 ,---,x p - 1 z 1 } 


(79) 


V 


■= {x°z\x°z 2 ,---,x°z p - 1 } 


(80) 



and 



for Vi, v 2 , ■ ■ ■ , v p and v , respectively. 

Proof. The proof of Proposition 15 is straightforward: It follows from (|62l) and (|63l) . For 
Proposition 16, we need to pass from u(p) to its subalgebra su(p). A basis for the Lie 
algebra su{p) of SU{p) is provided with the set {X a Z b : a,b = 0, 1, ■ ■ • ,p-l}\{X°Z }. 
Then, in order to prove Proposition 16, it suffices to verify that the p + 1 sets (or classes) 
V , Vi, • ■ -, V p -u and Z constitute a partition of {X a Z b : a, b = 0, 1, ■ ■ • ,p-l}\{X°Z } 
and that the p — 1 operators in each set commute one with each other. Proposition 16 
takes its origin in a remark RSI according to which the rank of su(p) is p — 1 so that the 
case of p + 1 sets containing p — 1 commuting operators occurs as a limiting case. The 
decomposition (1731) , also valid for s/(p, C), was first derived in [9] in connection with the 
determination of finest gradings of Lie algebras of type A p _i. It is little known that a 
decomposition of type (1731) was conjectured almost three decades ago [|59l for the more 
general case where p is replaced by a prime power (see also [60]). □ 
Example 3. For the purpose of clarifying the production process of the sets Vi (for i = 
0, 1, ■ ■ • , p), let us consider the case p = 7 ■$$ j = 3). Equations (|74l) -(l80l) give 



V 

V-i 

V 2 

v 3 

Vi 

v 5 

v 6 

V 7 



{(01) 


(02) 


(03) 


(04) 


(05), ( 


{(10) 


(20) 


(30) 


(40) 


(50), ( 


{(11) 


(22) 


(33) 


(44) 


(55), ( 


{(12) 


(24) 


(36) 


(41) 


(53), ( 


{(13) 


(26) 


(32) 


(45) 


(51), ( 


{(14) 


(21) 


(35) 


(42) 


(56), ( 


{(15) 


(23) 


(31) 


(46) 


(54), ( 


{(16) 


(25) 


(34) 


(43) 


(52), ( 



(06)} 
(60)} 
(66)} 
(65)} 
(64)} 
(63)} 
(62)} 
(61)} 



(81) 
(82) 
(83) 
(84) 
(85) 
(86) 
(87) 
(88) 
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where (ab) is used as an abbreviation of X a Z b . 

At this stage, it should be stressed that decompositions of type (I73H80I ) are especially 
useful for the construction of mutually unbiased bases [|40i |56l . Along this vein, the 
common eigenvectors of each of the p + 1 subalgebras v ,v 1 ,- ■ ■ ,v p give rise to p + 1 
mutually unbiased bases. Unfortunately, finding a general formula for the Lie brackets of 
each pair of the Cartan subalgebras is a difficult problem for which we have no answer. 
Counterexample 1. For d = 4 -^ j = 3/2 (=>- a, b = 0, 1, 2, 3), Proposition 15 is valid 
but Proposition 16 does not apply. Indeed, the 16 unitary operators u ab corresponding to 

ab = 01, 02, 03, 10, 20, 30, 11, 22, 33, 12, 13, 21, 23, 31, 32, 00 (89) 

are linearly independent and span the Lie algebra of U(A) but they give only 3 disjoint 
sets, viz., {(01), (02), (03)}, {(10), (20), (30)} and {(11), (22), (33)}, containing each 3 
commuting operators, where here again (ab) stands for X a Z b . However, it is not possible 
to partition the set (l89l) in order to get a decomposition similar to (|73T) . Nevertheless, 
it is possible to find another basis of u(A) which can be partioned in a way yielding a 
decompostion similar to (1731) . This can be achieved by working with tensorial products 
of the matrices X a Z b corresponding to p — 2. In this respect, let us consider the product 
u ai b! ® u a2b , 2 , where u aibi with i = 1, 2 are Pauli operators for p = 2. Then, by using the 
abbreviation (aibia 2 b 2 ) for u aibl ® u a2 b 2 or X ai Z bl ® X a2 Z b2 , it can be checked that the 
5 disjoint sets 

{(1011), (1101), (0110)} (90) 

{(1001), (0111), (1110)} (91) 

{(1010), (1000), (0010)} (92) 

{(1111), (1100), (0011)} (93) 

{(0101), (0100), (0001)} (94) 

consist each of 3 commuting unitary operators and that the Lie algebra sw(4) is spanned 
by the union of the 5 sets. It is to be emphasized that the 15 operators (I90H94I ) are underlaid 
by the geometry of the generalized quadrangle of order 2 [30]. In this geometry, the five 
sets given by (I90H94I ) correspond to a spread of this quadrangle, i.e., to a set of 5 pairwise 
skew lines [|30l . 

The considerations of Counterexample 1 can be generalized in the case d := did 2 ■ ■ ■ d e . 
e being an integer greater or equal to 2. Let us define 

uab ■= u aibl ® u a2b2 <g> • • • <g> u aebe A := a 1 ,a 2 ,---,a e B := b u b 2 , ■ ■ ■ , b e (95) 
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where u aibl ,u a2 b 2 , ■ ■ ■ , u ae b e are generalized Pauli operators corresponding to the dimen- 
sions d±, d 2 , ■ ■ ■ , d e respectively. (The operators uab are elements of the group P dl ® 
Pd 2 ® • • • <8> Pd e ■ We follow AH by calling the operators uab generalized Dirac operators 
since the ordinary Dirac operators correspond to P 2 (g> P 2 .) hi addition, let g 1; g 2 , • • • , g e 
be the g-factor associated with di, d 2 ,- ■■ ,d e respectively (g, := exp(27ri/d,-)). Then, 
Propositions 13, 14 and 15 can be generalized as follows. 
Proposition 17. The operators uab are unitary and satisfy the orthogonality relation 

T^£{did2-d e ) \u ab ua'B'J = did 2 ■ ■ ■ d e 5 a ,a> 5 b ,b> (96) 

where 

5 A, A' '■= S ai ,a'^a2,a' 2 " " " $a e ,a' e $B,B' '■= ^fei^&a,^ ' ' ' $b e ,b' e - (97) 

The commutator [uab,ua'b']- and the anti-commutator [uab, Ua'B']+ of uab and ua'B' 
are given by 

[uab, u A 'b']t = I II Qj "' =F II Qj ^ J ] u A "B" (98) 

with 

A" := oi + oi, a 2 + a' 2 , • • • , a e + a! e B" := b x + b[, b 2 + b' 2 ,---,b e + b' e . (99) 

The set {u A b : A,B E TLd^ ® ^d 2 ® ' ■ • ® ^d e } of the d\d\ ■ • ■ d\ unitary operators uab 
form a basis for the Lie algebra u(did 2 ■ ■ ■ d e ) of the group U{d\d 2 ■ ■ ■ d e ). In the special 
case where d\ — d 2 — ■ • ■ — d e — p with p a prime integer (or equivalently d = p e ), we 
have [uab, ua>b'}- = if and only if 

e 

Y^ cijb'j - bja'j = (mod p) (100) 

3=1 

and [uab, ua'b']+ = if and only if 

6 1 

2_. a jb'j ~ bjafj = -p (modp) (101) 

3=1 

so that there are vanishing anti-commutators only if p = 2. For d = p e , there exists a 
decomposition of the set {uab ■ A, B E Z® 6 } \ {/} that corresponds to a decomposition 
of the Lie algebra su(p e ) into p e + 1 abelian subalgebras of dimension p e — 1. 
Proof. The proof of (|96l) - (1101l) is based on repeated application of Proposition 13. For 
d = p e , we know from H24l [25l [33j that the set {u AB ■ A, B E Zf e } \ {1} (consisting of 
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p 2e — 1 unitary operators that are pairwise orthogonal), which provides a basis for su(p e ), 
can be partioned into p e + 1 disjoint classes containing each p e — 1 commuting operators. 
Therefore, there exists a decompostion of su(p e ) into a direct sum of p e + 1 subalgebras of 
dimension p e — 1. (There is a one-to-one correspondence between thep e + 1 subalgebras 
and the p e + 1 mutually unbiased bases in C p< \) □ 

5 Closing remarks 

Starting from an abstract definition of the Heisenberg-Weyl group, combined with a polar 
decompostion of SU(2) arising from angular momentum theory, we have analysed in 
a detailed way the interelationship between Weyl pairs, generalized Pauli operators and 
generalized Pauli group. The interest of these developments for the unitary group U(d), 
d arbitrary, have been underlined with a special emphasis for a decomposition of su(d) 
when d is the power of a prime. We would like to close with two remarks. 

In arbitrary dimension d, the number of mutually unbiased bases in C d is less or equal 
to d+ 1 [|24ll56l . Proposition 17 suggests the following remark. To prove that the number 
of mutually unbiased bases in C d is d + 1 for d arbitrary amounts to prove that it is 
possible to find a decomposition of the Lie algebra su(d) into the direct sum of d + 1 
abelian subalgebras of dimension d — 1. Therefore, if such a decomposition cannot be 
found, it would result that the number of mutually unbiased bases in C d is less than d + 1 
when d is not a prime power (cf. Conjectures 5.4 and 5.5 by Boykin et al Il60l0 . 

The Pauli group or discrete Heisenberg-Weyl group Pd = HW(Za) plays an impor- 
tant role in deriving mutually unbiased bases in finite-dimensional Hilbert spaces. We 
know that the concept of mutually unbiased bases also exists in infinite dimension IJoTfl . 
In this connection, the infinite or ordinary Heisenberg-Weyl group HW(M) might be of 
interest for constructing mutually unbiased bases in infinite-dimensional Hilbert spaces. 
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